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We have discovered a new antiferromagnetic phase in TmN^E^C by neutron diffraction. The 
ordering vector is Qa = (0.48, 0, 0) and the phase appears above a critical in-plane magnetic field of 
0.9 T. The field was applied in order to test the assumption that the zero-field magnetic structure at 
Qf = (0.094, 0.094, 0) would change into a c-axis ferromagnet if superconductivity were destroyed. 
We present theoretical calculations which show that two effects are important: A suppression of 
the ferromagnetic component of the RKKY exchange interaction in the superconducting phase, 
and a reduction of the superconducting condensation energy due to the periodic modulation of the 
moments at the wave vector Qa- 

PACS numbers: 74.70.Dd, 75.25.+Z, 74.20.Fg 
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The interplay between magnetism and superconductiv- 
ity is inherently of great interest since the two phenom- 
ena represent ordered states which are mutually exclusive 
in most systems. Therefore, the borocarbide intermetal- 
lic quaternaries with stoichiometry (RE)Ni2B2C have at- 
tracted great attention since the publication of their dis- 
covery in 1994 as they exhibit coexistence of mag- 
netism and superconductivity if the rare earth (RE) is ei- 
ther Dy, Ho, Er or Tm. The magnetic moments in these 
compounds are due to the localized 4/ electrons of the 
rare-earth ions. The 4/ and the itinerant electrons are 
coupled weakly by the exchange interaction, resulting in 
the indirect Ruderman-Kittel-Kasuya-Yoshida (RKKY) 
interaction between the 4/-moments. Thus the RKKY 
interaction, which is decisive for the cooperative behav- 
ior of the magnetic electrons, depends on the state of the 
metallic ones. 

The borocarbides have a tetragonal crystal structure 
with space group (14/mmm) ||, and TmNi2B2C has a 
superconducting critical temperature T c = 11 K and a 
Neel temperature T N = 1.5 K j5J. The crystalline electric 
field aligns the thulium moments along the c axis, and the 
magnetic structure has a short fundamental ordering vec- 
tor Q F = (0.094,0.094,0) with several higher-order odd 
harmonics |^,^|. In the magnetic structures detected in 
any of the other systems, the moments of the rare-earth 
ions are confined to the basal plane, and they have short 
wavelength antiferromagnetically ordered states. For ex- 
ample, they are commensurate with a propagation vec- 
tor Q = (0,0,1) for RE = Ho and Dy, or incommen- 
surate with Q w (0.55,0,0) for RE = Gd, Tb, Ho and 
Er J5j. An especially tight coupling between magnetism 
and superconductivity has been clearly demonstrated in 



TmNi2B2C, where a magnetic field applied along the c 
axis induced concurrent changes of the magnetic struc- 
ture and the symmetry of the flux line lattice JtJ . Similar 
effects have not yet been observed in any of the other 
borocarbides. 

One question that immediately arises, and which 
we believe is important for the general understanding 
of the interaction between superconductivity and mag- 
netism, is why the long-period magnetic ordering found 
in TmNi2B2C is stable. Bandstructure calculations on 
the normal state of non-magnetic LUNI2B2C predict a 
maximum in the conduction-electron susceptibility x(q) 
at q « (0.6,0,0) supported by Fermi-surface nesting 
H[| . The RKKY interaction is proportional to the mag- 
netic susceptibility of the electron gas, and the position 
of its maximum, determining the magnetic ordering vec- 
tor, is expected to be nearly the same as in x(q). This is 
in agreement with the experimental findings for many of 
the magnetic borocarbides, but not for TmNi2B2C. 

The use of the band structure of the normal state ig- 
nores that the magnetic interactions are mediated by a 
superconducting medium. The BCS ground state is a 
singlet and apart from the interband scattering the elec- 
tronic susceptibility is therefore zero at q — and zero 
temperature. It will recover its normal-state value only 
when q > 10 where £ is the coherence length of the 
superconductor, as first pointed out by Anderson and 
Suhl @. This raises the possibility that a local maxi- 
mum in the susceptibility is shifted from zero in the nor- 
mal phase to a non-zero value of q in the superconducting 
phase. In the case of the free electron gas, this maximum 
is very shallow and occurs at a relatively large q. If, how- 
ever, interband effects (or umklapp processes in the free 
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electron model) are included the situation is changed. We 
shall write the RKKY interaction in the superconduct- 
ing phase as J(q) = I[xo(q) + X« (<?)]■ Xo(o) is the con- 
tribution from the intraband scattering near the Fermi 
surface, which is sensitive to the superconducting energy 
gap. We have determined this function numerically in the 
zero temperature limit by using a linear expansion in q of 
the band electron energies near the Fermi surface. A sim- 
ple fit to the result is XoO) = 0.99 q [q + 1.5^ _1 ]- x Xq (0), 
where £ = nA/(Hvp) is the coherence length. The su- 
perconducting energy gap is negligible compared to the 
band splittings, and the contribution Xu(q) from the in- 
terband scattering is therefore not affected by the state 
of the conduction electrons. At small q we may assume: 
Xuil) — Xo (0)( Q! — Aq 2 ), where A is a constant normal- 
state quantity, which is expected to be some few times 
(27r/a) -2 , where a is the lattice parameter. The local 
maximum in J(q), at q = in the normal phase, will in 
the superconducting phase appear at qo « [(4/3)j4£] -1 ' 3 . 
This is the same dependence on £ as found by Anderson 
and Suhl, but the coefficient does not depend explicitly 
on kp and is somewhat smaller. With A(2tt/o) 2 lying 
between 1 and 10, qo assumes a value between 0.084 and 
0.18 times 2n/a, which is consistent with the magnitude 
of the observed ordering vector, Qp = 0.13(27r/a). 

The existence of a local maximum in the RKKY inter- 
action is not sufficient for stabilizing a magnetic structure 
with that wavelength. The total free energy of the entire 
system, including the condensation energy of the super- 
conducting electrons, has to be minimal. The presence of 
a magnetic modulation at Q may influence and weaken 
the superconductivity, for instance through a reduction 
of the density of states at the Fermi surface due to the su- 
perzone energy gaps created by the periodic modulation 
of the moments. This effect depends strongly on Q and 
may be an important factor in the competition between 
alternative magnetic structures. 

In this letter we report the results of neutron diffrac- 
tion studies on TmNi2B2C augmented by theoretical cal- 
culations of the competition between the magnetic and 
superconducting states. With the objective of suppress- 
ing superconductivity while the c-axis components of the 
magnetic moments are still ordered, a magnetic field was 
applied in the basal plane. This allows a study of the 
magnetic state as it would be in the absence of super- 
conductivity. The experimental results are described in 
detail below. Briefly we found that the system does 
not form a ferromagnetic state when superconductiv- 
ity is suppressed, as one might expect from the argu- 
ments presented above, but instead it relaxes into an in- 
commensurate antiferromagnetic state with Q = Qa — 
(0.482, 0, 0), resembling the magnetic ordering most com- 
monly found in the borocarbides. 

The experiments were performed on a 2 x 2 x 0.2 mm 3 
single crystal of TmNi2B2C grown by a high-temperature 
flux method and isotopically enriched with n B to en- 



hance the neutron transmission |ll|,|l2j. The sample was 
mounted on a copper rod in a dilution refrigerator insert 
for the low-temperature measurements at T < 1.7 K, and 
on a standard sample stick for measurements above 1.7 
K. In both cases the sample was oriented with the a- 
and &-axes in the scattering plane and placed in a 1.8 T 
horizontal-field magnet aligned along the a axis. Mea- 
surements were performed with applied fields up to the 
maximum of 1.8 T, and at temperatures between 20 mK 
and 6 K. The neutron diffraction experiments were per- 
formed at the TAS7 triple-axis spectrometer on the cold 
neutron beam line at the DR3 research reactor at Ris0 
National Laboratory. The measurements were carried 
out with 12.75 meV neutrons, pyrolythic graphite (004) 
monochromator and (002) analyzer crystals, and in an 
open geometry without collimation. The effective beam 
divergence before and after the sample was 60 and 120 
arc minutes, respectively. 
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Fig. 1 Normalized integrated intensities versus temperature of the 
field-induced magnetic peaks at 1.2 T (circles), 1.4 T (triangles) 
and 1.8 T (squares). The linear fit to the 1.8 T data shows the 
determination of T N (B). Inset: Scan along the [hOO] direction at 
T = 100 mK and 1.8 T, showing the field-induced satellite peaks at 
Q A = (0.482, 0, 0) around the (0 0) and (2 0) nuclear reflections. 
The peak intensity of the (2 0) reflection is 800. 

In the inset to Fig. 1 is shown the result of a scan 
along [10 0] at T = 100 mK and B = 1.8 T. In ad- 
dition to the nuclear (2 0) Bragg reflection, the figure 
shows satellite peaks at (0 0) + Qa and (2 0) — Qa 
with Qa = (0.482,0,0). These satellite peaks are not 
observed in zero applied field neither below nor above 
T N , and they do not appear until the field is larger than 
0.9 T at 100 mK. Additional magnetic satellites were de- 
tected at (1 1 0) ± Q A and (0 2 0) ± Q A No higher order 
harmonics of the field-induced satellites were observed, 
and their intensities are consistent with the magnitude 
and direction of the magnetic moments remaining equal 
to 3.8 /j, b and being essentially parallel to the c axis. At 
all temperatures and fields, where the Qa peaks were ob- 
served, they stayed resolution limited. Likewise no field 
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or temperature dependence of Qa was found. Quite re- 
markably, the field-induced peaks are only observed for 
the wave vector Qa || B and not for Qa -L B, hence the 
lowering of the in-plane four-fold symmetry of the sys- 
tem produced by the field has a direct consequence on 
the direction of Qa- 
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Fig. 2 Experimental and theoretical phase diagram of TmNi2B2C 
in a magnetic field along [10 0]. The medium-gray area denotes 
the region where both the Qa and the Qp reflections are present. 
In the dark-gray area only the Qa reflections were observed, up 
to the maximum applied field of 1.8 T. The light-gray area de- 
notes the region where the long tail of low-intensity magnetic scat- 
tering at Qa is still observed. The squares denote the measured 
phase boundary between the Qa phase and the paramagnetic one, 
T N (B), determined by the procedure described in the text. The 
open circles denote the upper critical field determined by trans- 
port measurements [ta|. The solid lines are the theoretical phase 
boundaries. The dashed line is the calculated Neel temperature 
of the Qa phase had the metal stayed in the normal state. The 
thin line labeled B® 2 (T) is the estimated upper critical field if the 
magnetic subsystem is neglected. 

Concurrently with the appearance of the field-induced 
peaks the intensity of the zero-field magnetic reflections 
with scattering vector Qp decreases and finally vanishes 
at 1.4 T and 100 mK. Between 0.9 and 1.4 T the magnetic 
structures at Qf and Qa coexist. The length of Qp 
does not change for applied magnetic fields up to 0.9 
T. Above 0.9 T a small reduction of \Qf\ is observed, 
simultaneously with the appearance of the field-induced 
magnetic reflection at Qa- The reduction is at the most 
3%, just before the peaks vanish at 1.4 T. 

In the main body of Fig. 1 we show the temperature de- 
pendence of the integrated intensity of the field-induced 
magnetic reflections for three different values of the ap- 
plied fields. This shows that the intensity of the peaks 
increases with increasing field at a constant temperature. 
The results at the different values of the field show quali- 
tatively the same temperature dependence, a rapid linear 
decrease when the temperature is above ~ 0.5 K and a 



cross-over into a long tail at about 2 K. At the maximum 
field of 1.8 T the tail extends up to a temperature of four 
times the zero-field value of T N . For each value of the 
field T N (B) is defined as the extrapolation to zero of the 
linear part of the integrated intensity, as shown in the 
figure for the case of 1.8 T. The experimental data are 
summarized in the phase diagram in Fig. 2. 

The theoretical phase diagram is calculated using the 
following parameters: 1) The crystal- field Hamiltonian 
of the Tm ions determined from experiments |l4|,|l5|]. 2) 
A phenomenological RKKY interaction J(q) with two 
parameters J{Qf) and J{Qa)- The normal-state value 
of J(0) w J{Q F )- 3) Abrikosov's formula JT^J for the 
condensation energy of the superconducting state 

F s - F„ = -{B° c2 (T) - BJ 2 [1.16 • 8^(2 K 2 - l)]" 1 , (1) 

where Bi is the internal magnetic field, corrected for the 
uniform magnetization and demagnetization effects. In 
the present system Bi is about 10% larger than the ap- 
plied field. B® 2 (T) is the upper critical field if the cou- 
pling to the magnetic electrons is neglected. Its depen- 
dence on T/T c is assumed to be the same as observed in 
the non-magnetic Lu borocarbide IlT]]. B® 2 (T — 0) has 
been used as a fitting parameter, and its final value of 6.5 
T (corresponding to £ = 71 A) is close to that obtained 
when scaling the values of B® 2 (0) m the non-magnetic Lu 
and Y borocar bides with T c , The other fitting parameter 
in Eq. (1) is the (renormalized) value of k which is found 
to be 6.3, close to that determined experimentally by Cho 
et al. (ll]] . 4) The coupling between the magnetic system 
and the superconducting electrons is described by two 
parameters: a) The Anderson-Suhl reduction of ^7(0) in 
the superconducting phase at zero temperature, which is 
found to be close to the value of J7"(0) itself. The ther- 
mally excited quasiparticles implies that the reduction 
is smaller at finite temperatures. This effect is included 
in the calculations, b) The suppression of the supercon- 
ductivity, which is assumed to be due to the superzone 
energy gaps near the Fermi surface produced by the mag- 
netic ordering at Qa- The density of states at the Fermi 
surface is reduced proportionally to the sizes of the en- 
ergy gaps, which by themselves are proportional to the 
amplitude of the magnetic modulation psf . The reduc- 
tion of the density of states causes a decrease of B® c2 (T) 
in Eq. (1), which is estimated to be at most 40%. 

The model calculations account reasonably well for the 
experimental phase boundaries in Fig. 2, and equally well 
for the phase diagram when the field is applied in the 
c direction, as measured by Eskildsen et al. M. The 
large anisotropy between the experimental upper critical 
field along the a and along the c axis is not an inter- 
nal property of the superconducting electron system, but 
is due to the large difference between the magnetic a- 
and c-axis susceptibility of the Tm ions. Although the 
transformation of the Qf state into the c-axis ferromag- 
net did not occur, it is clear that the Anderson-Suhl 
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mechanism plays an important role for the behavior of 
TmNi2B2C, see also the discussion by Kulic et al. [ p^[ . 
The loss in the superconducting condensation energy at 
the field-induced transition to the normal state is com- 
pensated for by the gain in the RKKY exchange energy 
deriving from the uniform magnetization because of the 
sudden increase of J7(0). This mechanism explains the 
large reduction of the upper critical field shown in Fig. 
2, and the even larger reduction in the c-axis phase di- 
agram. Most significantly, it explains why the Qf state 
is stable up to a field of 1 T applied along the c axis, i.e. 
as long as the system stays superconducting. In case of a 
normal behavior of J(q), the c-axis field would destroy 
this phase almost immediately. 

Against intuition, the short wavelength magnetic mod- 
ulation in the Qa phase affects the superconductivity 
state much more strongly than the long wavelength Qf- 
modulation, but it should be noticed that the wave- 
lengths of the two structures are both much shorter than 
£. The model predicts that in a normal magnetic system 
the Qa phase would be the stable one at low temper- 
atures and fields. However, the superzone energy gaps 
produced by the Qyi-modulation disturb the nesting fea- 
tures of the Fermi surface close to this wave vector, caus- 
ing a strong suppression of the superconductivity in the 
Tm system. This effect is also found in the measure- 
ments on e.g. Ho based borocarbides, where supercon- 
ductivity is suppressed when the magnetic system enters 
the Q — (0.55, 0, 0) phase and is regained when the sys- 
tem leaves this phase at a yet lower temperature p0| . Er 
borocarbide is the only one where the (0.55, 0, 0)-phase 
does not seem to strain the superconducting properties 
severely, although the upper critical field p2| shows sim- 
ilarities with the observations in the Tm system. The 
reason that the Qf phase does not appear in, for in- 
stance, the Er system may simply be that the RKKY 
interaction, being proportional to [(g — 1) J] 2 , is stronger 
in Er than in Tm making the energy difference between 
the two different magnetic phases too large in comparison 
with the superconducting condensation energy. 

Improvement of the present theory is planned in order 
to account for the inhomogeneitics in the superconduct- 
ing order parameter due to the flux lines. It is difficult to 
understand how the low-intensity long tail of the Qa re- 
flections may extend up to a temperature which is twice 
the estimated transition temperature of the QA-phase in 
a normal system (the dashed line in Fig. 2). The phase 
may survive within the normal core of the flux lines, but 
not very far above the normal-phase transition tempera- 
ture. Further experiments are planned at higher fields in 
order to test the theoretical predictions, e.g. that T N (B) 
coincides with the maximum in the in-plane upper critical 
field, to investigate in more detail the magnetic scattering 
in the long-tail regime, and to search for the Qyi-phase 
in the c-axis phase diagram. 

In this letter we have proposed that the long wave- 



length magnetic structure of TmNi2B2C below 1.5 K 
owes its very existence to the fact that the material is 
superconducting. This was shown experimentally by ap- 
plying an in-plane magnetic field, which in this Ising-like 
system primarily affects the superconducting state. We 
have found that the magnetic system enters a new phase 
at a critical field of 0.9 T where the order is at short 
wavelength, Qa = (0.48,0,0). We have presented the- 
oretical calculations showing that the interplay between 
superconductivity and magnetism in the borocarbides is 
governed by two mechanisms: 1) a suppression of the 
ferromagnetic component of the RKKY exchange inter- 
action in the superconducting phase, and 2) a reduction 
of the superconducting condensation energy from the pe- 
riodic modulation of the moments at the wave vector Qa- 
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